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Abstract. 

In 2001, using a large spectroscopic dataset from an extensive monitoring campaign, we discovered that the peculiar Of star 
HD 108 displayed extreme line variations. This strange behaviour could be attributed to a variety of models, and an investigation 
of the high energy properties of HD 108 was needed to test the predictions from these models. Our dedicated XMM-Newton 
observation of HD 108 shows that its spectrum is well represented by a two temperature thermal plasma model with kTi ~ 
0.2 keV and kT2 ~ 1.4 keV. In addition, we find that the star does not display any significant short-term changes during the 
XMM-Newton exposure. Compared to previous Einstein and ^05Ar detections, it also appears that HD 108 does not present 
long-term flux variations either. While the line variations continue to modify HD 108's spectrum in the optical domain, the X- 
ray emission of the star appears thus surprisingly stable: no simple model is for the moment able to explain such an unexpected 
behaviour. 

Thanks to its high sensitivity, the XMM-Newton observatory has also enabled the serendipitous discovery of 57 new X-ray 
sources in the field of HD 108. Their properties are also discussed in this paper. 
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1. Introduction 

Some O-type stars are still challenging astronomers many 
decades after their discovery. HD 108 is a good example of such 
an object. In the past, this star has been the target of several in- 
vestigations, but apparently, none could arrive at a consensus 
on the exact nature of this peculiar star. Consequently, various 
theories have successively been proposed: HD 108 has been 
classified as a short-term binary (having even survived a super- 
nova event according to some authors); it has been proposed to 
be a single star experiencing wind variability and/or harbour- 
ing a disc and jets; and finally it was suggested to be a long- 
term binary (see Naze et al. 2001, hereafter Paper I, and ref- 
erences therein). In Paper I, we presented a 30 year campaign 
of optical spectroscopy dedicated to HD 108. We showed that 
the behaviour of HD 108 was not that of a classical short- or 
long-term SB 1 binary. The star also did not present any short- 
term variations. But our extensive campaign clearly indicated 
the peculiar characteristics of this star: tremendous line varia- 
tions on the timescale of decades. The hydrogen and He i lines 
change from strong P Cygni profiles to simple absorptions, 
and a few other emission lines apparently follow the same be- 
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haviour Comparing with all data available, we noted that these 
variations were recurrent, on a timescale of a few decades. Such 
a continuous decline of H i and He i lines was recently discov- 
ered by Walborn et al. (2003) in another Of?p star, HD 191612. 
The line variations of this star are strikingly similar to HD 108, 
but occur on a shorter timescale. In this context, we note that 
HD 191612 has now returned to a high emission state: a spec- 
trum of this star covering the Ha Une, taken in October 2003 at 
the Observatoire de Haute-Provence, is almost identical to that 
of August 1997 presented by Walborn et al. ( 2003 ) . 

The nature of these Of?p stars is still unknown. The most 
popular models to explain their peculiar behaviour involve the 
presence of a compact companion on an elliptical orbit, or re- 
quire the star to be a single rapid rotator (see e.g. the review 



in Walborn et al. I2003t . It was also suggested that these stars 
could be transition objects, explaining their small number (only 
3 are known in our Galaxy). In this context. X-ray observations 
represent an important opportunity to better understand these 
objects, since stars reveal at high energies the most exotic pro- 
cesses taking place in their vicinity. X-ray data are especially 
well suited to test the possibility of accretion processes linked 
to the presence of a compact object or a colliding wind inter- 
action in a binary. We thus decided to observe HD 108 with 
the XMM-Newton satellite, in the hope that its high sensitiv- 
ity could provide us with definitive answers on the nature of 
HD 108. 

In this paper, we describe in Sect. 2 the observations used 
in this study. The optical and X-ray data of HD 108 will then 
be successively analysed in Sect. 3. The remarkable sensitivity 
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Fig. 1. EPIC hardness ratio map of the whole field around 
HD 108. Three energy bands were used to create this color im- 
age: red corresponds to 0.4-1.0 keV, green to 1.0-2.0 keV and 
blue to 2.0-10.0 keV. HD 108 is the bright source at the center 
of the field. 



of XMM-Newton also enabled the serendipitous discovery of 
many fainter sources during our observation. These sources de- 
tected in the field of HD 108, their possible counterparts, their 
hardness ratios (HRs), their variability, and their spectral char- 
acteristics will be discussed in Sect. 4. Finally, we will con- 
clude in Sect. 5. 



2. Observations 

2.1. X-ray data 

HD 108 was observed with XMM-Newton in the framework of 
the guaranteed time of the Optical Monitor Consortium during 
revolution 494, on Aug. 21 2002, for approximately 35 ks. The 
two EPIC MOS cameras (Turner et al. 2001 1 were used in full 
frame mode, and the EPIC pn instrument (Striider et al. 200]} 
was operated in extended full frame mode. A thick filter was 
added to reject optical light. 

We used the Science Analysis System (SAS) software ver- 
sion 5.4.1. to reduce the EPIC data. These data were first pro- 
cessed through the pipeline chains, and then filtered. For EPIC 
MOS, only events with a pattern between and 12 and passing 
through the #XMMEA_EM filter were considered. For EPIC 
pn, we kept events with flag=0 and a pattern between and 
4. To search for contamination by low energy protons, we ex- 
amined the light curve at high energies (Pulse Invariant channel 
number > 10000, ~E>10keV, and withpattern=0). A very short 
flare was detected and we thus discarded the time intervals with 
a high energy count rate larger than 0.2 cts s ' (for EPIC MOS) 
and 0.64 cts s ' (for EPIC pn). The resulting useful exposure 
times are 35.7 ks and 28.9 ks for EPIC MOS and pn, respec- 
tively. Further analysis was performed using the SAS v5.4.1. 
and the FTOOLS tasks. The spectra were analysed and fitted 
within XSPECvl 1.0.1. 

The first two figures show images of the field from the com- 
bined EPIC instruments. Figure ^ presents a three colour im- 
age, in which numerous hard X-ray point sources are clearly 
seen. Figure Ogives the identification of the point sources de- 
tected in the field (see Sect. 4 and TableO. 

The RGS data were processed through the 'rgsproc' 
pipeline SAS meta-task. This routine determines the location 
of the source spectrum on the detector (see e.g. den Herder 
et al.^01), and performs the extraction of the source spec- 
trum as well as of the background spectrum from a spatially 
offset background region. Finally, the task computes tailored 
response matrices for the first and second order RGS spectra. 



2.2. Additional optical data 

In 2001, 2002 and 2003, we continued our extensive spec- 
troscopic survey of HD 108 with the Aurelie spectrograph at- 
tached to the 1.52 m telescope of the Observatoire de Haute- 
Provence (OHP). The detector was a 2048x1024 CCD EEV 
42-20#3, whose pixel size is (13.5fj.mf-. The exact wavelength 
ranges and mean S/N are given in Table ^ together with the 
dispersion. All the data were reduced in the standard way us- 
ing the MIDAS software developed at ESO. The spectra were 
normalized by fitting splines through carefully chosen contin- 
uum windows, the same as the ones used in Paper I. 



Table 1. Summary of our CCD spectroscopic observations. 
indicates the total number of spectra obtained during a run. 



Date 


Wav. range 


N 


S/N 


Disp. 


R 










(A mm"') 




Sep. 2001 


3680-4130 


2 


160 


16 


7000 


Sep. 2001 


4080-4540 


2 


160 


16 


7000 


Sep. 2001 


4450-4900 


6 


200 


16 


7000 


Sep. 2001 


6350-6760 


2 


360 


16 


11000 


Sep. 2002 


4450-4905 


2 


250 


16 


7000 


Sep. 2002 


6350-6760 


1 


260 


16 


11000 


Oct. 2003 


4450-4915 


5 


40-350 


16 


7000 


Oct. 2003 


6350-6750 


1 


300 


16 


11000 



3. HD108 

3.1. Latest optical data 

After a reduction consistent with that of the data from previous 
years (see Paper I), we found that HD 108 had still not reached 
its minimum emission stage. The H i. He i. Si m, and C m lines 
continued to decline. The average equivalent widths' (EWs) 
were 1.66+0.02Ain 2001 for the Hyline; 1.02+0.03Ain2001, 
1.11+O.OlA in 2002 and 1.22+0.04A in 2003 for He i /14471; 
0.72+O.OlA in 2001, 0.72+0.02A in 2002 and 0.72+0.04A in 
2003 for He ii /14542; 2.01+0.03A in 2002 and 2.08+0. 12A 
in 2003 for H/3. We show in Fig. the evolution of the EWs 
with time: the changes are striking. We also note that the Hj3 
line may now also present variations on a shorter timescale: 
in 2001 and 2003, the absorption apparently increased in the 
last spectra of the runs. The spectral type of the star has also 
continued to change through the years, as we had predicted in 
Paper I: it was 08Ifpe in 2001-2002 and 08.5fpe in 2003. 

In 1997, 2001, 2002 and 2003, we also obtained spectra of 
HD 108 in a region containing the Ha line. This line too under- 
went dramatic variations, as is shown in Fig.|4] its EW changes 
from -5.34A in 1997 to -1.37A in 2003! The profile is now 
asymetric, with a small absorption in the blue wing. We can use 
this Ho- EW to estimate the mass-loss rate of the star (Lamers & 
Leitherer ll993l . But first, we have to correct the observed EW 

' In Table 5 of Paper I, the EW of H/3 should be read -1.25 A in 
1987 and 1.49±0.07A in 2000. 
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HD108 : all EPIC in 0.4-10 keV 
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Fig. 2. EPIC image of HD 108 in the range 0.4-10 keV. The detected sources are labelled. The image has been binned by a factor 
of 50, to obtain a pixel size of 2.5". 



for photospheric Hi and He ii absorptions. To this aim, we de- 
cided to use the observed Ha EW of 081 and OS. 51 stars listed 
in Conti (19741. However, the only value available in Conti's 
paper, 2.91 A, is that of HD 167771, which is currently known 
to be a binary star. Since this could affect the determination of 
the EW of its Ha line, we rather use the average EW as mea- 
sured by Conti ( 1974.) on the spectra of single 07-9.51 stars 
(HD 34656, HD 163800, HD 16429, HD 218915). This value 
is 3. 06 A. Using the method outlined in Lamers & Leitherer 
| |1993j with the physical parameters appropriate to HD 108 (see 
Sect. 3.2.1), an observed V magnitude of 7.4 and color (B- 
V) of 0.18 mag (see Paper 1), we then find that the mass-loss 
rate log(M) varies from -5.4 to -5.5 between 1997 and 2003, 
whatever the photospheric EW correction chosen. These val- 
ues of M are in between the extreme estimations from the lit- 



erature (see below). Note however that the method of Lamers 



& Leitherer J1993> assumes that the stellar wind is spherically 
symmetric and homogeneous. The former hypothesis might not 
be valid in the case of HD 108 (see discussion below). 

In 2001, we observed HD 108 in the violet range between 
3680 and 4130 A. No O iii /1/13755,3760 emission was detected. 
These lines would be seen in emission if the N iii /1/14634,4640 
emissions were produced through a Bowen mechanism in- 
volving the O m and N m transitions at 374A. But since these 
O III lines are on the contrary in absorption, this mechanism is 
clearly not at work in the wind of HD 108. The Niii emission 
lines are then most probably explained by a combination of the 



processes described in Mihalas J1973> . i.e. excitation of certain 
autoionising states of N iii followed by dielectronic recombina- 
tions (Mihalas . 1973 J and pumping of adequate N iii transitions 
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Fig. 4. Evolution of the red spectrum of HD 108 with time. In the 1997 data, the Ha line was in pure emission while in the 21st 
century data, it exhibits a clear P Cygni profile, whose emission component decreases with time. 
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Fig. 3. Average EWs for all observing runs: He ii /1 4200 is 
represented by filled triangles, Hy by distorted squares. He i 
A 447 1 by crosses. He ii A 4542 by filled circles, and H/3 by open 
triangles. The errors (evaluated from the data dispersion) are 
very small, roughly the size of the symbols used in the figure. 

by the intense continuum radiation of the star in an expanding 
atmosphere (Swings . 1948j . 

3.2. X-ray data analysis 

EPIC spectra and lightcurves of HD 108 were extracted over a 
circular region centered on the star and of radius 50" for MOS 
and 37'.'5 for pn (because of the presence of a gap nearby). 



Since a fainter source is present near HD 108, we did not 
choose an annular background region, but rather a nearby cir- 
cle devoid of sources. The EPIC lightcurves did not present any 
significant variation. The spectra of each instrument were fit- 
ted separately in XSPEC, but as they gave the same properties, 
within the errors, we decided to fit simultaneously all EPIC 
data available. We fixed the interstellar absorbing column to 
A^h(ISM) = 3.4 X 10^1 cm-2 (Diplas & Savage lT99^ . The best 
fit was then obtained with a sum of two absorbed optically thin 
plasma models (mekal, Kaastra ri992> of temperatures kT - 
0.2 and 1.4 keV (see Fig.|5}- A fit by a differential emission 
measure model (c6pmekl. Sing et al. ll996> does not improve 
the quality of the fit but confirms the existence of two domi- 
nant temperatures. We note the presence, at ~6.6 keV, of the Fe 
K-Hne. 



Figure|6] shows the first order RGS 1 and RGS2 spectra of 
HD 108. Unfortunately, the signal to noise ratio of individual 
lines is not sufficient to perform a detailed quantitative analysis 
of the RGS line spectrum of HD 108. However, the emission 
lines that are seen in these spectra clearly indicate that the low 
energy (<2 keV) X-ray emission arises in a thermal plasma. 
We identified the strongest emission lines by comparison with 
the line list of the SPEX plasma code (Kaastra et al. 2002 1. The 
temperatures of maximum emissivity of the lines seen in the 
RGS spectra span a range from ~ 2 x lO*" to ~ 10^ K. These 
values are compatible with the temperatures obtained with the 
spectral fits of the EPIC data. We therefore included the RGS 
data in a simultaneous fit of all the X-ray data of HD 108 and 
the results are presented in Table |2] Unlike the cases of f Pup 
(Kahn et al. l2Uon and 9 Sgr (Rauw et al. l2002a> . almost no 
lines are seen above ~ 20 A, where the spectrum is heavily 
absorbed by circumstellar and interstellar material. 
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Table 2. Spectral properties of HD 108. The spectra were fitted simultaneously on the data from the three EPIC cameras and the 
two RGS gratings. The fluxes are in the 0.4-10.0 keV energy range and the emission measures are given for a distance of 2.51 
kpc (Gies 1987 1. The models were of the type wabs * {absi * mekali + abs2 * mekalj), with the first absorbing column fixed to 
the interstellar value, A^//(ISM) - 3.4 x 10^' cm"^ and the unabsorbed fluxes listed below were dereddened only for this column. 
The absorption components absi and absi are either from neutral gas {wabs) or an ionised wind (see text). The quoted errors 
correspond to the 90% confidence intervals. 



Abs. 


Column 1 
10^2 cm-2 


keV 


cm"' 


Column2 
10^2 cm-2 


keV 


EM2 

cm 


Flux 

10-'- 


Unabs. Flux 

-2 -1 

erg cm s 


x'(d.o.f.) 


wabs 
wind 




20"-^- 


9.4x10^* 
3.7x10^* 


1 SO'""* 


1 Q7I.44 
1 54 'fi^ 


6.9x10^' 
5.8x10^5 


6.85 
7.00 


15.7 
17.6 


1.18 (852) 
1.15 (852) 



3.2.1. Wind modelling 

The material in the stellar winds of early-type stars is ionised 
by the stellar radiation. Several authors (e.g. Krolik & Kallman 
119841 Waldron ll984t have demonstrated that a comprehensive 
study of the X-ray throughput from massive stars requires a 
detailed modelling of the opacity from such an ionised wind. 
We have thus attempted to model the wind opacity of HD 108. 

We have considered the 10 most abundant elements (H, 
He, C, N, O, Ne, Mg, Si, S, Fe) and fixed their abundances 
to the solar ones (Anders & Grevesse 1989). In our model, 
the coUisional excitation, the photoionisation and the radia- 
tive and dielectronic recombinations determine the ionisation 
of the elements. CoUisional excitation rates were taken from 
Voronov (.1997.) . and the photoionisation cross sections from 
Verner et al. J1996> for the outer shells and Verner & Yakovlev 
([1995^ for the inner shells. Since innershell ionisation can lead 
to the ejection of several electrons, we used the Auger yields 
from Kaastra & Mewe {1993} to model this effect for the el- 
ements considered. Radiative and dielectronic recombination 
rates were taken from Verner & Ferland ( 1996 1 and Shull & 
Van Steenberg (I1982> . taking into account the corrections from 
Arnaud & Rothenflug ('1985) for the latter reference. 

We determined the ionisation structure of the stellar wind 
up to ~300 stellar radii (with a logarithmic spatial bin from 
5 R, to 321 /?,), using the approximations of Waldron ( 1984i 
for the radiation field (his eq. 12) and the opacities (his eqs. 
13 & 14). However, we chose to fix the velocity law and the 
wind temperature to v(r) = v^, (1 - 0.99^)"-^ and T„i„d(r) - 
To. + (To - Too)(^y '^ (with To = T,ff, Too =0.4*7^0, see 
Lamers & Morris 1994), respectively. The stellar parameters 
needed in these equations were taken from several sources: 
Howarth & Prinja ([1989) and Peppel ([1984) give = 17 Rq, 
M ~ 65 Mo and Te/f ~ 40 kK, while Howarth et al. ((1997) 
quote Voo = 1960 kms In addition, note that the stellar flux 
was taken from Kurucz's library of spectra for the above effec- 
tive temperature and log(§)=4.5 (the closest value available) 
and that we modelled the local emission as a simple blackbody 

There is no consensus in the litterature about the mass-loss 
rate of HD 108: values from 2xlO-^M0yr ' (Hutchings & Von 
RudlofflW) to 5xlO"^M0yr-' (FeiTari-Toniolo et al. rMTb 
can be found. On the other hand, using the Ha line EW as mea- 



HD 108 (EPIC MOS+pn) : wind model 



0.5 1 2 5 10 

channel energy (keV) 



Fig. 5. The EPIC MOS and pn spectra of HD 108, with the best 
fit model with wind absorption (see Table |2) superimposed. 
Black lines are for EPIC MOSl data, fit, and residuals; light 
gi-ey is for EPIC MOS2 and dark grey for EPIC pn. 

sured on our spectra from 1997-2002, we found an intermedi- 
ate value (see above). To test the effects of changing mass-loss 
rates on the wind ionisation, we decided to run our model for 
the two extreme mass-loss rates, and in the rest of the paper, 
we will refer to these models as 'wind-low' and 'wind-high', 
respectively. 

We show in Fig.0the resulting wind optical depth for both 
models at five stellar radii, compared to the interstellar one. 
If we scale the optical depth t at any radius r in the wind by 
Nwmd(r) - ^ p{r')dr\ we find that the resulting opacity is 
very similar whatever the radius or the mass-loss rate consid- 
ered in the model. We thus implemented a unique wind opacity 
model in XSPEC (we chose the wind-high opacity at 5 R,) and 
then tried to fit the XMM-Newton spectra using this absorb- 
ing model. The results are presented in Table |5] We see that 
within the errors, the absorbing columns are very similar com- 
pared to the case of neutral gas. The temperatures are slightly 
larger and the Emission Measures (EMs) lower, compared to 
the case of the neutral gas, but the;^f^ value is just slightly lower. 



6 



Naze et al.: HD 108: the mystery deepens with XMM-A'ewton observations 



0.025 

^ 0.02 

''^ 0.015 
w 

o 

o 0.01 



^ 0.005 
o 

S5 




10 15 

Wavelength (A) 

Fig. 6. The RGSl and RGS2 spectra of HD 108 between 6 
and 23 A, along with the best fit wind model of Table |2l The 
strongest lines are identified. 



For an ionised wind, the two fitted columns coiTespond respec- 
tively to positions in the wind at 86 and 37 /?» for the largest 
M (5xlO"''MQyr"'), or 7 and 3 i.e. relatively close to the 
star, for the intermediate one (M =3.5xlO"^Moyr"', a value 
appropriate for the epoch of the XMM-Newton observation if 
the wind is spherically symmetric, as assumed in our model). 
Note that such high columns could not be reached with a low 
value of the mass-loss rate like M =2x lO^^Moyr"' . Such a low 
mass-loss rate actually produces a much lower optical depth 
than the interstellar matter (see Fig.Q. The absorption by such 
a wind would thus not lead to a significant signature at soft X- 
ray energies, contrary to what is observed: very low mass-loss 
rates are not compatible with the X-ray data. 

3.3. Comparison with previous X-ray observations 

The XMM-Newton observatory is actually not the first to detect 
HD 108 in X-rays. The star had already been observed twice in 
the last 25 years. In January 1979, during a 1800 s exposure, the 
Einstein satellite detected an X-ray source, 2E0003. 4-1-6324, 
at the position of HD 108. The IPC count rate for this source 
was 2.8+0.5 10"^ cts s ' (2E catalog, Harris et al. [T994t or 
3.2+0.6 10^2 cts s^' (Chlebowski et al. 1989^. Twelve years 
later, /?(9&4r reobserved this region during 558 s in the course 
of the All-Sky Survey. This time, the X-ray counterpart of 
HD 108 was named IRXS J000605.0+634039 (Voges et al. 
I1999j and presented a count rate of 6.1 + 1.3 10"^ cts s"'. 

Using the spectral properties of Table |2] we can predict 
the count rates expected with the IPC and PSPC_C instruments 
if these properties have not changed since then. Whatever the 
absorption model, the XMM-Newton data convert into an IPC 
count rate of 2.6±0.5 10"^ cts s"' and a PSPC.C count rate of 
4.8+2.9 10"^ cts s"'. The observed count rates are thus compat- 
ible, within the eiTors, to what could be expected on the basis of 
the XMM-Newton data: we observe at most a marginal decrease 
of the flux. Contrary to the optical domain, there is thus appar- 




-4 



3 4 6 8 10 

Energy (in keV) 



Fig. 7. Logarithm of the optical depth in the X-ray domain from 
5 Rt outwards. The solid line corresponds to M =2xlO"^M0 
yr"', and the dotted one to M = 5xlO"^Mo yr"'. For compar- 
ison, we also show with a dashed line the optical depth of the 
neutral hydrogen column due to the interstellar medium (ISM): 
Nh(ISM) = 3.4 x 10^1 cm-- (Diplas & Savage . 1994.) . 



ently no dramatic change at the X-ray energies on a timescale 
of decades. 



3.4. Tine nature ofHD 108 

Before we consider the various models for HD 108, let us first 
highlight the most important results of our XMM-Newton in- 
vestigation. The XMM-Newton observations of HD 108 have 
revealed a thermal spectrum, of characteristic temperatures 0.2 
and 1 .4 keV. Moreover, the X-ray flux of HD 108 appears rather 
stable on both short timescales (of the order of the observation 
length) and longer ones (of the order of decades). 

To better understand the nature of the system, we can first 
compare the X-ray luminosity of HD 108, corrected for the in- 
terstellar absorption, to its bolometric luminosity. In the ROSAT 

band (0.1-2.0 keV), we find log(-^^) ~ -6.15, whereas we 
would expect a value of -6.80 from the relations of Berghofer 
et al. C1997V HD 108 is thus slightly overluminous in X-rays, 
with a factor of 4 between the expected and observed lumi- 
nosities. But this difference is not as large as e.g. in some 
coUiding-wind binaries, and our observed ratio is in fact just 
at the upper limit of the dispersion of the data of Berghofer 
et al. L^^^VLbol relations were also studied more recently by 
Moffat et al. (2002) in the 0.5-10 keV energy range, using the 
Chandra data of NGC 3603: HD 108 is then again amongst the 
stars at the upper limit of the dispersion. 

The emission measure (EM) of the entire wind (assuming 
spherical symmetry) is EM - J ngnndV ~ p^P'dr. 
For the intermediate mass-loss rate M =3.5x IO^^Mq yr ' and 
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with p{r) = j;^y we find EM ^3.9 X W cm^l By compar- 
ing this value to the fitted EM?, (see Table|2}, we can derive the 
volume filling factor of the hot gas. It is approximately 10"^ for 
the lower temperature component and ~ 2x10 for the higher 
temperature one. These filling factors are in agreement with 
values derived by Kudritzki et al. (1996 1 or Feldmeier et al. 
\\991t for other early-type stars. However, we may note that 
for the largest M considered above, the factor would be 200 
times lower and become, then, more difficult to reconcile with 
current X-ray emission theories and observations. 

3.4.1 . A binary system with a compact object 

In Paper I, we had proposed that HD 108 could be an eccentric 
binary consisting of the Of star and a compact object with an 
orbital period of several decades. In such a system, the compact 
object could accrete matter from its companion near periastron. 
At that moment, if the X-ray luminosity produced by wind ac- 
cretion becomes sufficient, the ionisation structure of the O-star 
wind could be severely altered, and e.g. elements like hydrogen 
and helium could become completely ionised. The wind emis- 
sions of H I and He i would thus disappear, thereby explaining 
the observed optical line variations. In such a model, we would 
expect in the XMM-Newton observation (1) a large X-ray lumi- 
nosity capable of affecting the ionisation of the wind and (2) 
large variations of the X-ray flux, on orbital timescales, corre- 
lated with the modulation of the optical emissions. However, 
the XMM-Newton observation revealed a completely different 
behaviour. The ISM-corrected X-ray luminosity measured by 
XMM-Newton is rather modest (L™°** ~ lO^"* erg s"' in the 0.4- 
10 keV energy range) and appeared stable compared to previ- 
ous X-ray observations made in the last 20 years. It is true that 
the ionisation of the wind can alter its radiative line driving, 
thereby modifying the mass-loss rate and thus probably the X- 
ray emission, but a constant X-ray luminosity could then only 
be explained by a delicate equilibrium between the changing 
wind absorption and the varying emission measure, which is 
unlikely. 

Let us consider what would be required to produce a change 
in the ionisation of the stellar wind as suggested by the opti- 
cal data. Tarter et al. ('1969 1 found that the ionisation equilib- 
rium in a stellar wind illuminated by an X-ray source depends 
mainly on the parameter ^ = where Lx is the X-ray lumi- 

"'~x 

nosity, n the local number density, and rx the distance to the 
X-ray source. Hatchett & Mc Cray (1977i futher introduced 
a dimensionless parameter q - ^'-^j^, where nx is the num- 
ber density at the position of the X-ray source and D the or- 
bital separation between the primary star and the X-ray source. 
Interestingly, q is essentially set by the velocity law of the stel- 
lar wind. The changes in the intensity of H i and He i emissions 
in the optical spectrum of HD 108 suggest that the properties 
of the wind are changing over a large volume. Therefore, the 
results of Hatchett & McCray ( 119771 their Figs. 1 & 2) indicate 
that q should be about ~ 1. Considering a number of differ- 
ent situations, Kallman & McCray (I1982> found that log(^) at 
the transition between neutral and ionised H and He should be 
in the range 1 .0 to 1 .5 depending on the wind density and the 



spectral properties of the X-ray source. These considerations 
suggest that ^ = ~ 10 - 30. 

On the other hand, the X-ray luminosity due to wind ac- 
cretion onto a neutron star can be estimated from Davidson 
& Ostriker (1973 ) as a function of the orbital separation D. 
Although we have shown that the X-ray spectrum of HD 108 
can be described by an optically thin thermal plasma model, 
and this is not the kind of spectral shape expected for an 
accretion-powered X-ray emission, let us assume for the mo- 
ment that the intrinsic Lx of HD 108 could entirely be assigned 
to the accretion by the compact companion. After deredden- 
ing the spectra by both interstellar and circumstellar absorp- 
tions, we get Lx ~ 0.75 - 1.5 10^"* erg s"' in the 0.4-10. keV 
energy range. These values correspond to orbital separations^ 
of 6 - 8/?. for M =3.5x 10~''Mo yr"', and this results in 

£ _ 4a-Lyl.3m/jv<f) ^ Q 16 — 32 
q M 

The intrinsic X-ray luminosity of HD 108 as determined 
with XMM-Newton is thus probably not sufficient to signifi- 
cantly alter the ionisation of the stellar wind and hence trigger 
the changes seen in the optical spectrum. We note however that 
a low X-ray luminosity does not completely rule out the pres- 
ence of a compact companion, since the wind accretion pro- 
cess can be inhibited under certain circumstances (Stella et al. 
I1986> . but that, if present, the physical conditions are such that 
this companion cannot explain the optical variations. 

3.4.2. A confined wind model 

A second type of model proposed is a rapid rotator with a shell 
and/or disc and/or jets (e.g. Walborn et al. 2003 1. If a disc 
is present in the system, it can only be seen nearly face-on: 
no double-peaked profiles have ever been reported for HD 108 
and spectropolarimetric studies also favor a face-on geometry. 
Since there is no modulation of the optical data on timescales of 
a few days that could be associated with the rotation of the star, 
the rotation and disc axes need to be nearly aligned. In such a 
system, long-term optical variations could result from preces- 
sion effects e.g. in a wide binary, but no change in the width of 
the lines, and hence of v sin /, is detectable in our spectra. 

The formation of a confined outflow could be linked to 
the presence of a strong magnetic field confining the wind to- 
wards the equatorial regions of the star. Such a magnetically 
confined wind could provide a straightforward explanation for 
the hard X-ray component seen in the EPIC spectra of HD 108. 
Models of such confined winds predict a strong shock at the 
interface between the winds that are deviated from both hemi- 
spheres and the equatorial disc (see e.g. Babel & Montmerle 
.1997. ud-Doula 2002 1. A change in the overall mass-loss rate 
could easily account for the changes of the intensity of optical 
recombination emission lines such as Ha but would also im- 
ply a change of the total wind emission measure proportional 
to AM^ which should then lead to a significant correlated vari- 
ation of the X-ray luminosity that is not observed. 



- These orbital separations correspond to orbital periods of 
16-25 days if the orbit was circular (although this is unlikely to be 
the case for HD 108). 
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3.4.3. Other models and summary 

Another model associates HD 108 with transition objects 
(Paper I, Walborn et al. 120031 and references therein) such as 
Ofpe/WN9 stars or Luminous Blue Variable stars (LBVs). No 
particular characteristics in the XMM-Newton data argue in fa- 
vor of such a model, but these data do not provide arguments 
against either However, as in the previous model, it is difficult 
to understand why a tremendous change in the wind structure 
where emission lines arise is not correlated with a similar dra- 
matic variation at the wind position where the X-ray emission 
takes place. 

Finally, we note that the second plasma temperature 
{kT2 - 1.4 keV) is rather high. Most single O-type stars 
display thermal X-ray emission with kT < \ keV and higher 
temperature emission is often attributed to a colliding wind 
phenomenon in a binary system. A colliding wind scenario in a 
wide (eccentric) O+O binary could also account for the slight 
X-ray overluminosity of HD 108. However, it is not clear how 
the changes in the optical spectrum could be related to this 
scenario. Also, we emphasize that there is no obvious evidence 
for such a binary system in the spectrum of HD 108. 

All models are now challenged by the XMM-Newton ob- 
servations. Their context does not permit us to understand how 
the large variations seen in the optical domain could lead to a 
stable X-ray luminosity. This apparent dichotomy could well 
be explained by a difference in the formation region of the X- 
ray emission and the H i/He i emissions, e.g. X-rays could arise 
in a disc while the emission lines could be formed in a jet, 
but it is difficult to reconcile this with the fact that a change 
in any physical parameter of the star would affect only one of 
the two features. For example, an increase in the mass-loss rate 
could only lead to a stable X-ray luminosity if the resulting 
additional X-ray absorption and emission of the wind compen- 
sate exactly. Such a delicate balance between emission and ab- 
sorption was invoked by Owocki & Cohen ( 1999) to explain 
the empirical Lx/Lbol relation of O-type stars. This formalism 
holds for a mass-loss driven by radiation pressure. However, 
since no large changes of the magnitude of HD 108 have been 
reported, it seems unlikely that a change in the mass-loss rate 
of the star could be triggered by a change in the overall lumi- 
nosity and it seems therefore excluded that the mechanism of 
Owocki & Cohen could explain the lack of X-ray variability of 
HD 108. More than ever, HD 108 is a puzzle for astronomers. 

4. X-ray sources in the field 

4.1. List of sources 

The X-ray data of HD 108 also reveal a lot of additional point 
sources in the field surrounding the O star, as can be appreci- 
ated in Fig. n To study these discrete X-ray sources, we ap- 
pUed the detection metatask edetectjchain simultaneously to 
the data from the three EPIC cameras. We used three energy 
bands: 8=0.4-1.0 keV, M=1.0-2.0 keV and H=2.0-10.0 keV. 
We eliminated the false detections mainly due to background 
fluctuations by rejecting the sources with a likelihood < 5 in 
any detector and/or a combined likehood < 30. This procedure 



left 58 point sources in the field. They are shown in Fig.|2and 
their properties are presented in Table|3] by order of increasing 
RA. 

Their background-subtracted and vignetting-corrected 
count rates in each of the EPIC cameras are also presented in 
this table. Note that the SAS detection tasks do not apply a cor- 
rection for out-of-time events or dead times. For the EPIC pn 
data used in extended full frame mode, this implies an over- 
estimate of the exposure time (and thus an underestimation 
of the count rates) by 2.3%. Since this amount is well within 
the errors, we choose not to correct the pn count-rates pro- 
vided by edetect -chain. The sensitivity limit is ~0.6 cts s"' for 
EPIC MOS data and -1.6 cts s ' for EPIC pn data. The hard- 
ness ratios for the pn data are also indicated in Table 13 The 
hardness ratios are defined as HR\ - (M - S)/(M + S) and 
HR2 ^ (H - M)/(H + M), where S, M, and H correspond to 
the count rates in the 0.4-1.0, 1.0-2.0 and 2.0-10.0 keV band, 
respectively. As was already visible from Fig. ^ most of the 
sources are very hard, probably due to the large absorbing col- 
umn in the Galactic plane where HD 108 lies {In - 117.93°, 
bu = +1.25°). 

4.2. Identification 

We have compared our source list to the ROSAT detections in 
the same field. Four catalogs are available: the WGA catalog 
(White et al.|2)00 PSPC) and three catalogs from the ROSAT 
teams (IRXH and 2RXP /?C>5Ar consortium 2000 1 and 200jb; 
IRXS, Voges et al. 1999. & .2000). Only one source, IRXS 
J000605. 0+634039, was previously known in the field. It is at 
13'.'7 from XMMU J000603.2+634046 (#37) and at 13'.'3 from 
HD 108. Although this distance is slightly larger than the po- 
sitional error of the ROSAT source, 11", there is little doubt 
that the ROSAT source corresponds to HD 108, since there is 
no other bright XMM-Newton source in the region. 

To search for optical counterparts, we have cross-correlated 
our source list with the Simbad catalog, the 2MASS All-Sky 
Data Release, the USNO Bl.O and the GSC 2.2 catalogs. To 
determine an optimal radius of correlation, we adopted the 
method of Jeffries et al. (I1997> as used in Rauw et al. ( 20 02bt . 
i.e. fitting the cumulative distribution of the number of detected 
sources as a function of the correlation radius by expression (1) 
of Jeffries et al. ( I1997> . The optimal correlation radius^, i.e. the 
radius that includes the maximum of true correlations and the 
minimum of spurious correlations, is found to be around 2". 
We list the counterparts to the X-ray sources found this way 
in Tabled together with the separation between the X-ray and 
the visible sources, and some basic information about the visi- 
ble sources. 

Source XMMU J000428.5+633121 (# 7) is actually ex- 
tended, as can be seen from comparison with the shape of the 
neighbouring sources (see Figs. ^ and |2li. This source is not 



Strictly speaking, a unique correlation radius cannot be used with 
XMM-Newton data, since the point spread function is degrading at 
large off-axis angles. The radius found here should be considered as a 
weighted average. 
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associated with any GSC/2MASS/USNO counterpart: with its 
rather large HR, it may be an extragalactic source. 

The source XMMU J000643.6+633912 (# 48) corresponds 
to BD+62°2368, an F7V star If we compare the observed col- 
ors of this star (with B, V listed in Simbad and 7, H, K from 
the 2MASS catalog) to typical colors of an F7V star (Schmidt- 
Kaleril982 Kenyon & Hartmann 1995 Bessell & Brett 1988 1, 
we can estimate its reddening and its probable distance. Using 
the reddening law from Cardelli et al. ( 1989) with =3.1, we 
found an ^(B - V) in the range 0.02-0.07, and a distance of 
181-193 pc. 

The source XMMU J000623.9+633927 (#40) also corre- 
lates with a bright star, BD+62°2365, but no spectral type has 
been quoted for this star in the literature so far However, its 
2MASS colors and a visible spectrum obtained at OHP in Oct. 
2003 apparently favor a G2III spectral type for this star. In a 
similar manner as above, we derived for this object a redden- 
ing ^(B - V) ~ 0.18 and a distance of 540 pc. 

On the other hand, several sources listed in Table|3]present 
rather high HRs. This result does not appear completely 
surprising in the light of the log(A^) - log(5) relations of 
Mushotzky et al. ('2000) and Giacconi et al. ( 1200 H . since these 
relations predict that a large fraction of the detected sources 
should correspond to extragalactic objects even accounting for 
a total Galactic column density of ~0.5xlO^^ cm"^ along the 
line of sight'*. Finally, there are a few known early-type stars in 
the field in addition to HD 108 and EM* CDS 1: BD+62°2362 
(B7III), ALS 6022 (B9III), ALS 6028 (B), EM* CDS 3 (B), 
BD-h62°2369 (A5V) and BD+63°2105 (A9II). None of them 
is detected in X-rays, but such late B and A stars are usually 
not bright X-ray emitters. We may also note the non-detection 
in X-rays of the variable star S V*S VS 1455 and of the relatively 
bright ALS 6011 (F5I). 

4.3. Variability 

We have analysed the lightcurves of the 8 brightest point 
sources (i.e. with at least 50 cts in EPIC MOS and 100 cts in 
EPIC pn in the 0.4- 10.0 keV range). The count rates in each bin 
were background-subtracted using annuli around the sources or 
close-by circles. The eflFective time bin lengths were calculated 
by taking into account the 'good time intervals' defined in Sect. 
2. The resulting lightcurves were analysed by Kolmogorov- 
Smirnov and tests. We also ran a modified probability of 
variability test (Sana et al. l2003> on these sources' event lists. 
In the HD 108 field, no source was found significantly variable 
over the duration of our observation. 

4.4. Spectral properties 

We also extracted the spectra of these brightest sources. We 
generated response matrix files (rmf) and ancillary response 
files (arf) using the SAS tasks rmfgen and arfgen. The spec- 

We estimated this value from the extinction map of Schlegel et al. 
<1998> . These authors caution that their maps may not be accurate for 
Galactic latitudes below \b\ ^ +5°, but they can nevertheless provide a 
first estimate of the total extinction. 



tra were then binned to reach a minimum of 10 cts per channel. 
Finally, we analysed the background-corrected spectra within 
XSPEC. Due to strong noise at very low and very high en- 
ergies, we have discarded energy bins below 0.4 keV and 
above 10 keV. We have fitted the spectra with an absorbed 
mekal model or an absorbed power-law. For each source, 
we have fitted separately EPIC MOSIh-2 and EPIC pn data, 
but as they gave similar results - within the errors - we fi- 
nally fitted all three instruments simultaneously and we list 
the parameters of the best-fit models in Table |4] Note that 
XMMU J000643.6H-633912 (# 48), which corresponds to the 
F7V star BD+62°2368, is only well fitted by a mekal model, 
as could be expected for a coronal source. For the aver- 
age distance of 187 pc estimated above, the X-ray luminos- 
ity of this source is 8.2x10^^ erg s"\ which corresponds to 

log(-^^) = -5.1. A mekal model also provides a good fit 
for XMMU J000623.9H-633927 (# 40), and for the distance 
of 540 pc estimated above, its X-ray luminosity is 4.5x10^'' 

erg s"', and that results in log(^^^) = -5.6. On the contrary, 
the large absorbing column, the large HRs, and the good fit 
by a power law of XMMU J000445. 1+633550 (#12), XMMU 
J000457.6+634211 (# 14) and XMMU J000701. 7+634151 (# 
54) suggest that these 3 sources may actually be background 
objects. The case of XMMU J00051 1.7+634018 (# 19) and 
XMMU J000537. 1+634027 (# 23) is different: they are also 
better fitted by a power law, but their absorbing column is small 
and their fluxes seem rather low for close X-ray binaries. 

5. Conclusion 

We have obtained an XMM-Newton observation of HD 108 and 
its surroundings. The peculiar star HD 108 was found to present 
a two-temperature spectrum, and did not show any significant 
short-term variations during the exposure. These observations 
are also compatible with a stable X-ray emission since 1979. 
In parallel, we have continued our extensive optical monitoring 
of the star, and discovered that HD 108 continues to present 
dramatic line variations in the optical domain. The lack of sig- 
nificant changes in the X-ray emission compared to the opti- 
cal data is a puzzle. In such a context, no simple model is for 
the moment capable of explaining the overall behaviour of one 
of the most peculiar O stars of the sky: any variation of the 
physical properties (magnetic field, mass-loss rate) of HD 108 
is expected to have an impact on the whole wind, not only on 
the formation region of the optical emission but also on that of 
the X-rays. Long-term monitoring of HD 108, preferentially in 
a multi-wavelength campaign, is thus necessary to eventually 
understand this star. 

In addition, 57 new X-ray sources were also discovered in 
the field of HD 108, and we present here their characteristics 
(count rate, HRs, plus lightcurve analysis and spectral fits for 
the brightest ones). Only two correspond to rather bright stars, 
both of them being late-type stars with coronal X-ray emission. 
With large HRs and/or large absorbing columns and power law 
fits, several sources might be background objects. On the other 
hand, none of the B stars of the field was convincingly detected 
in X-rays. 



Table 3. Characteristics of sources in the field of HD 108. The hardness ratios refers to EPIC pn, and the counterparts are within 2" of the X-ray sources. For each catalog, 
the column 'Nr' shows the number of counterparts lying in this 2" region, and 'd' gives the distance in arcseconds between the X-ray and the optical source, if unique. The 
quoted error bars on the count rate represent ±lcr poissonian standard deviation. 
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36 


XMMU J000559.5+634149 


1.04+0.23 


0.99+0.24 


3.25+0.48 


0.53+0.21 


0.00+0.15 



























37 


XMMU J000603.2+634046'' 


99.9+1.8 


105.+2. 


333.+4. 


-0.10+0.01 


-0.50+0.01 


1 


0.7 




7.5 


1 


0.7 


6.9 


6.9 


6.8 


1 


0.7 


38 


XMMU J000606.4+633347 


1.17+0.29 


1.22+0.31 


2.12+0.50 


-0.65+0.19 


-1.00+1.34 



























39 


XMMU J00062 1.9+634321 


0.68+0.21 


0.57+0.20 


0.88+0.33 


1.00+0.73 


0.51+0.32 



























40 


XMMU J000623.9+633927'' 


2.06+0.28 


2.20+0.32 


9.51+0.70 


-0.47+0.06 


-0.76+0.15 


1 


1.1 


9.7 




1 


0.7 


8.3 


7.8 


7.7 


1 


0.8 


41 


XMMU J000624.8+633737 


1.54+0.29 


1.33+0.27 


3.77+0.53 


0.72+0.31 


0.51+0.12 



























42 


XMMU J000629.5+634106 


1.02+0.26 


1.30+0.27 


2.93+0.53 


0.05+0.19 


-0.25+0.23 


1 


0.5 


14.6 


15.5 


1 


0.4 


12.8 


12.3 


12.2 


1 


0.2 


43 


XMMU J000632.3+634127 


0.68+0.23 


0.53+0.18 


2.47+0.48 


0.50+0.24 


-0.17+0.21 



























44 


XMMU J000638.4+634902 


0.73+0.29 


0.67+0.27 


2.25+0.58 


-0.16+0.23 


-0.54+0.49 


1 


1.5 


14.2 


16.0 


1 


0.4 


13.8 


13.2 


13.0 


3 




45 


XMMU J000639.8+633627 


1.55+0.31 


1.07+0.28 


3.33+0.52 


0.68+0.25 


0.36+0.13 



























46 


XMMU J000640. 1+634255 


0.63+0.21 


0.74+0.21 


2.00+0.46 


0.32+0.28 


0.02+0.25 



























47 


XMMU J000642.9+633308 


0.67+0.33 


1.07+0.38 


3.39+0.69 


0.31+0.24 


0.14+0.21 



























48 


XMMU J000643.6+633912' 


2.95+0.31 


2.55+0.31 


13.5+0.1 


-0.76+0.04 


-1.00+0.07 


1 


1.8 




10.9 


1 


1.7 


9.3 


9.1 


9.1 


1 


1.8 


49 


XMMU J000644.2+633739 


1.48+0.30 




4.47+0.60 


0.63+0.17 


-0.03+0.14 



























50 


XMMU J000647.2+633602 


0.78+0.29 


0.63+0.25 


1.00+0.40 


1.00+0.55 


-0.17+0.38 



























51 


XMMU J000647.5+634947 


1.19+0.38 


1.20+0.35 


1.21+0.52 


0.97+0.31 


-0.52+0.51 



























52 


XMMU J000658.1+633221 


3.11+0.65 


2.85+0.53 


7.02+1.01 


-0.71+0.45 


0.97+0.04 



























53 


XMMU J000701.3+635130 


4.38+1.25 


5.12+0.73 

































54 


XMMU J000701 .7+634151 


7.99+0.62 


8.33+0.63 


24.0+1.3 


0.95+0.08 


0.71+0.03 



























55 


XMMU J000705.7+634447 


1.18+0.25 


1.66+0.33 

































56 


XMMU J000707.8+634721 


0.89+0.29 


1.03+0.34 


1.55+0.62 


0.27+0.42 


-0.04+0.47 



























57 


XMMU J000713.2+633439 


1.98+0.44 


1.39+0.40 


5.73+0.88 


0.55+0.22 


0.28+0.14 



























58 


XMMU J0OO718.8+633857 


1.69+0.37 


1.06+0.32 


3.01+0.64 


-0.66+0.17 


-0.44+0.69 


1 


1.5 


16.0 


18.3 


1 


1.6 


13.2 


12.5 


12.3 


1 


1.6 



Extended source 

The emission-line star EM* CDS 1 (B-type) is at 8'/5 from this sotirce. 
^ = HD 108 (the X-ray position is at ff.' 8 from ttie optical one given in Simbad). 
_ = BD+62°2365 (spectral type G2in, see text) 
^ ^ = BD+62°2368 (spectral type F7V) 
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Table 4. Spectral properties of the brightest sources of the HD 108 field-of-view. The spectra were fitted simultaneously on the 
data from the three EPIC cameras, except for # 19 and 23, where there exist only MOS data. The observed fluxes are in units 
10 erg cm"^ s ' and in the 0.4-10.0 keV energy range. Models with bad fit quality (x^ > 2) are not Usted. The quoted errors 
correspond to the 90% confidence intervals. 



Absorbed power-law model 

r Flux Ar^(d.o.f.) 
keV 



Nil 



Absorbed mekal model 

kT Flux ^^(d.o.f.) 



10^2 cm- 



Var. ? 



12 
14 
19 
23 
40 
48 
54 



1 27'-''* 
64" '"' 
0.02|j;"' 

"•"0.09 

61°''' 



1 63'-'^ 



3.43|- 
82j°o 
42J00 
23°-35 
39" '" 

0.32 
IQlOO 



1.06 
0.60 
1.08 
1.03 
0.13 
0.20 
1.72 



0.83 (74) 
0.91 (65) 
0.62 (47) 
1.04 (52) 
0.79 (51) 
1.47 (62) 
0.74 (95) 



0-84j:^^ 

0.13 

"•0. 

29" ''2 

1 201-^* 



1 Q92.64 
^•^^1.43 



2 06^'* 

^•""l.76 

7 472-82 

1 11 '■'•2 

1 40'-'*' 
10'" 



1 60''^ 



1.16 
0.58 
1.18 
1.01 
0.13 



0.83 (74) 
0.90 (65) 
0.60 (47) 
1.03 (52) 
0.96 (51) 



1.69 0.72 (95) 



N 

N 

? 

N 
N 

N 
N 
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